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FOREWARD 


This report describes the work performed between September 
1988 to December 1988 under Research Grant Number NAG-1-685, 
entitled, "Use of Energy Accountancy and Power Flow Techniques 
for Aircraft Noise Transmission". This is the fifth progress 
report under this research grant in which most of the work done 
thus far addresses the problem of the characterization of the 
vibrational power flow. During this latest phase of the work, a 
Ph.D. Thesis was completed with the title, "Vibrational Power 
Flow in Thick Connected Plates", the Abstract of which will be 
published in the Journal of the Acoustical Society of America in 
the Technical Notes and Research Briefs section. The work in 
this thesis was also partially funded by an ONR Fellowship which 
paid for the stipend of the graduate student. 

The author would like to acknowledge the graduate students 
who assisted in the work on this project, the Department of Ocean 
Engineering and most importantly the financial support by the 
Structural Acoustics Branch of the NASA Langley Research Center. 
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ABSTRACT 


This report investigates the problem of the vibration power 
flow through coupled plate substructures when one of the 
substructures is subjected to a distributed load. The approach 
used is based on the mobility power flow method. Two types of 
distributed loading are considered. The first part of the report 
deals with a distributed mechanical load which is independent of 
the response of the structure. Power flow expressions are 
presented for this type of loading. In the second part of the 
report, the analysis for distributed load which is influenced by 
the motion of the structure, as in the case of acoustic 
excitation is presented. Results for power flow when the plate 
substructures are subjected to acoustic excitation will be 
published in future reports. 
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1 . 


INTRODUCTION 


The previous work using the mobility power flow (MPF) 
approach has, in general, been restricted to point excited 
structures [1-5]. However in the case of excitation of an 
aircraft fuselage distributed loading on the surface of a panel 
can be as important as excitation from directly applied forces at 
defined locations on the structure. The term "mobility power 
flow" is being used here to distinguish this technique from other 
techniques on power flow. In this approach, the vibrational 
power flow is determined using mobility expressions. There are 
other methods which are similar to either finite element analysis 
[6] or statistical energy analysis where the results are 
frequency or spatial averaged [7,8]. In this report the 
formalism of the power flow expressions for two coupled plate 
substructures with a distributed load on one of the substructures 
is presented. 

Two types of loadings are considered. They are a distributed 
force load of arbitrary shape and a distributed load created by 
an obliquely incidence acoustic waves on the surface of one of 
the plate substructures. The difference between these two types 
of distributed loads is that in the mechanical loading case the 
applied loading is independent of the response of the structure 
while in the case of the excitation from an incident acoustic 
wave, the loading is dependent on the response of the structure 
due to the presence of the scattered acoustic waves. In this 
latter type of loading coupled expressions must be developed for 
the excitation and response. This, to some extent, can be viewed 
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as an introduction into the development of expressions for the 
power flow (acoustical and vibrational) , using a mobility power 
flow approach, when coupling exists between a structure and the 
surrounding medium as in the case of the interior of an aircraft 
fuselage. 


2 . DISTRIBUTED FORCE 

The development of power flow and mobility expressions for 
the distributed force loading, with arbitrary dependence on the x 
and y coordinates, is very similar to the work reported 
previously for point loading [2,4], The differences are that the 
presentation of the mobility functions can be generalized to 
apply for any loading condition. Using the configuration and 
axes definition as shown in figure 1, the modal input mobility 
function for the junction, and the transfer mobility function 
between the junction and any point on the plate surface, both 
functions being defined as the response per unit applied moment 
for every mode m along the junction, are given by [2,4]. 
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where T m is the mode m component of the edge moment, k^and k 2 are 
defined by: 


2 2 2 

k, - 2k + k ; k„ = k 

1 x y 2 y 



u m is the mode m component of the plate surface displacement, 
is the mode m component of the angular displacement at the edge 
of the plate and D* and ph are the plate flexural rigidity and 
surface density respectively. In the above equations it is 
assumed that the two plate substructures are identical. 

Note that in eguation (2) the transfer mobility is a function 
of coordinate y. This is required since the load is distributed 
and defined everywhere along the y direction on the plate 
surface. 

The above modal mobility expressions are for a system with an 
applied edge moment, which represents one of the configurations 
for which mobility expressions are derived [2], The second 
configuration to be considered is the one which represents the 
external loading to the global structure. Considering first the 
distributed load, this can be decomposed into its modal 
components as follows: 

F(x,y) = F q f(x)f(y) (4) 
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F mn ~ F m F n ( 7 ) 

This load decomposition will be used to define the mobility 
functions for the second plate configuration. Therefore, the 
modal mobility for the input load, Mi mn defined as the ratio of 
mode (m,n) component of the transverse velocity on the surface of 
the plate, to an applied distributed surface load with the same 
distribution as the eigenfunction for mode m,n. is given by: 
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and u mn is the mode m,n component of the plate surface 
displacement. The expression for the transfer mobility, M2im 
defined as the ratio of the mode m component of a plate edge 
angular velocity, to the mode m component of an applied 
distributed load on the surface of the plate is given by; 



( 10 ) 


The evaluation of the above mobility functions is obtained 
from consideration of the equations of motion for a plate 
structure. 

The expression for the input and transmitted power are the 
same as in the case of point loading [2] except that all mobility 
functions are modal expressions; 
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m and n are respectively the mode number for the x and y 
directions, that is m is along the junction and n perpendicular 
to the junction, and a and b are the dimensions of the plates 
(figure 1) . 

For a point load, the force is in this case described 
mathematically by; 


F(x,y) = F 0 6 (x-x 0 ) S( Y-y 0 ) 


(13) 


and 
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where x Q and y 0 are the x and y coordinates of the point of 
application of the load. If the load is applied at the center of 
the plate (x Q = a/2 and y D = b/2) , then only odd values for m an 
n are allowed. Substituting in equations (ll) and (12) , the 
results obtained for the input and transferred power, when the 
load is in the center of the source plate, are shown in figure 
(2) . Comparing these results to those obtained in reference 
[13], the two sets of results are identical. This verifies the 
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formulation of the power flow expressions in the form shown in 
equations (11) and (12) . 

For a uniform distributed load described by: 

F(x,y) = F 0 (15) 


f (y) = i.o 


(16) 


and 
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The solution in this case follows in the same way as for the 
point load case. The results for the power flow are shown in 
figure (3) . 

Comparing the results for the power flow obtained for the 
point load with the results obtained for the distributed load, 
the following similarities and differences can be observed. 
First, the modes excited by the two types of loading are 
identical. This is expected since both loading conditions are 
symmetrical about the center of the source plate. Second, the 
general level of the power flow for the distributed load case 
decreases with frequency, while that for the point load does not 
decrease. The reason for this can be mathematically described by 
the inverse dependency of the modal components of the distributed 
applied load on the mode number for the case of the distributed 
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load (equation (17)). Physically this implies that the higher 
frequency modal components of the load are suppressed. The 
results for the distributed load would be similar to those for 
excitation by normal incidence acoustic plane waves, if the 
scattered pressure component is neglected. 


3. EXCITATION BY INCIDENT ACOUSTIC WAVES 
In dealing with this form of distributed excitation, the 
scattering of the incident acoustic wave due to the response of 
the receiving structure must be included in the analysis [9]. 
The scattered acoustic field is a function of the response of the 
structure, and therefore, the loading on the structure consists 
of two components; one component is from the incident pressure 
being blocked by the surface and the second component is the 
influence of the scattered acoustic wave. The component 
associated with the blocked incident pressure is independent of 
the motion of the structure except for the phase differences 
between structure response and surface motion. The component 
associated with the scattered pressure is a direct function of 
the motion of the structure and therefore the loading becomes 
structure response dependent. Therefore to deal with this type 
of distributed loading some modifications to the basic 
formulation of the expressions for the vibrational power input 
and power flow are required. To demonstrate the differences in 
the formulation, the case of a point load which is a function of 
the response velocity at the point of application for arbitrary 
structure configuration is discussed first. 
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Consider the configuration shown in figure 4, where the input 
load F is a function of the response velocity. The governing 
equations of motion for this configuration are as follows [1]; 

V = F + F a M 12 18. 


V A = “ f a m 3 


19. 


= F a M 2 + F M 2 i 


20 . 


but 

F = A +BV 21. 

Where A and B are two constants and all the other terms have 
the same definition as in previous reports in this series [1]. 
Equations (19) and (20) can be used to solve for the junction 
force F a in terms of the input load F, and then substituted into 
equation (18) to solve for the velocity V. In the case where F 
is a function of V, the result for V must be substituted back 
into the expression for F and then used to solve for F a and V A . 
A lternatively V and F A are expressed in terms of the constants A 
and B and the mobility functions. That is, 
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The input and transmitted power, and respectively, are 
then given by; 
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Applying this approach to the case of the acoustic 
excitation, the expressions for the input and transfer modal 
mobility functions for an input edge moment (equations 1 and 2) 
are still applicable. For the input mobility of an acoustic 
excited plate structure, consider an incident acoustic wave at an 
angle of 4 > to the normal to the plate, where the trace of the 
acoustic wave on the plane of the plate makes an angle 9 to the 
x-axis (figure 5) . The trace of the incident pressure on the 
surface of the panel is given by; 


, . _ -ik sin 9 sin 6 y -ik sin 6 cos 9 x 

p(x,y) « P i e J r J e J Y 25. 

where k is the acoustic wavenumber, and p^ is the magnitude of 
the incident pressure. 

The component of the input load associated with the blocked 
incident pressure is equal to twice the incident pressure [9] at 
the surface of the plate. The scattered component can then be 
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formulated based on the in vacuo response of the plate and is 
controlled by the amplitude of this response. The response that 
must be considered in determining the scattered component is the 
response of the global structure that is a simultaneous solution 
must be obtained for the forcing and the response of the global 
structure . 

Let the surface velocity of the receiver plate be described 

by; 


v(x,y) 


v(y) q 
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The scattered pressure can be obtained by equating the plate 
surface velocity to the particle velocity of the acoustic medium 
at the interface. Using a spatial transform to represent the 
response of the whole plate, the scattered pressure component for 
mode q is given by; 
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k z - ± J (k 2 - a 2 - p 2 ) 28 . 

where k is the acoustic wavenumber. The positive or negative 
sign is selected to satisfy the far field conditions. The total 
modal pressure decomposed into the wavenumber p,. acting on the 
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plate surface is given by the summation of the blocked component 
and the scattered component. 

P(£)m = 2 Pi(/3) m + P s (£)m 29. 


where 


oo 


a 


p l < «« - 1 

- f OT1 -jkxsin^ cosd 

2? i e J sm 

mnx 

. a J 

dx 

-oo 

o 




e -jky sin<£ sin0 



30. 


and 


P 08) - 

s m 


V</3) 


q 


nm 

[(-l) m e Jaa - 1] 

q7r 

(<-l) q e‘ J “- 1] “'’o do 

a 


2 

a 

r ^ 

mn 

2 ' 


. a . 


2 

a 

r q?r ' 

2 ' 

k 7ra 

z 

K ~ J 

L 3 - 

L 3 . 
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The terms of the summation for which q ^ m represent 
intramodal coupling which is negligible for light fluid loading 
provided the modal density of the structure is low. Equations 
(30) and (31) define the terms A and B in matrix equation (22) . 
To solve for the velocity of the source substructure; 


v <«» - »<«,, M <« m + 


T m «12 <«. 


where 


M(/3) 


m 




1 


* 

D 


r.2 

/■ ** 
m7r 

2 ' 

/3 + 




, a . 



phw 

* 

D 


32. 


33. 


15 


which is obtained from the plate equations of motion, and Mi 2 (£)m 
is the spatial fourier transform of equation (2) . From 
continuity of motion at the junction. 
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The integrand of this last expression contains a V(/3) m term 
and therefore when substituting equation (35) in equation (32) to 
solve for V(/9) , the solution can only be obtained numerically. 
From the solution for V(£), the input and transmitted vibrational 
power calculated. To evaluate the input power it is not 
necessary to inverse Fourier transform the expression for V(^) m 
and P(/3) m since an integral over the structure surface is 
required. This would be equivalent, by Parceval's Theorem, to an 
integral over the entire /? space. For the transmitted power this 
is obtained from the relationship; 
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and then summed for all modes m. 
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A complete solution to these equations has not been achieved 
and therefore no results will be presented in this progress 
report. Results for the L-shaped plate structure, including a 
parametric evaluation and comparison to the results for the point 
force loading, will be presented in a following progress report. 

4. CONCLUSION 

The mobility power flow formalism for two coupled structures 
with excitation from distributed loading, including acoustic 
excitation, on one of the components of the coupled structure has 
been discussed. For a distributed loading which is independent 
of the motion of the structure, the mobility power flow formalism 
is relatively straight forward and very much similar to that of 
point loading. For the case of acoustic excitation, where the 
magnitude of the excitation is controlled by the response of the 
global structure, the formalism is much more complex. Because of 
the coupling of the response to the excitation, the only method 
available for obtaining a solution is through a numerical 
approach. This is similar to conclusions reached by other 
authors who investigated the problem of acoustic excitation 
where, unless the mode shapes of vibration of the global 
structure are a priori known, the solution is always obtained 
numerically. This does not in any way limit the application of 
the mobility power flow method to multiple coupled structures, 
since the complexities arise mainly for the source structure. 

Results for the analysis presented in this report as applied 
to an L-shaped plate will be reported in the future together with 
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the experimental evaluation of the results for both point loading 
and distributed loading. 


REFERENCES 


1. J.M. Cuschieri, "Power Flow As A Complement to Statistical 

Energy Analysis and Finite Element Analysis", ASME 
Publication NCA-Vol 3, 58-61, (1987). 

2. J.M. Cuschieri, "Extension of Vibrational Power Flow 
Techniques To Two-Dimensional Structures", NASA Contract 
Report 181710, September, 1988. 

3. J.M. Cuschieri, "Parametric and Experimental Analysis 

Using A Power Flow Approach", NASA Contract Report 181990, 
February, 1990. 

4. J.M. Cuschieri, "Power Flow Analysis Of Two Coupled Plates 
With Arbitrary Characteristics", NASA Contract Report 
182033, June, 1990. 

5. R.J. Pinnington, R.G. White, "Power Flow Through Machine 

Isolators to Resonant and Non-Resonant Beams", Journal of 
Sound and Vibration, 75, 179-197, 1981. 

6. D.J. Nefske, S.H. Sung, "Power Flow Finite Element 

Analysis of Dynamic Systems: Basic Theory and Application 

to Beams" ASME Statistical Energy Analysis Symposium 
proceedings, NCA - Vol . 3, 47-54, 1987. 

7. L.D. Pope, J.F. Wilby, "Band Limited Power Flow Into 

Enclosures", Journal of the Acoustical Society of America, 
vol. 62 (4), 906-911, 1977. 

8. L.D. Pope, J.F. Wilby, "Band Limited Power Flow Into 

Enclosures IT", Journal of the Acoustical Society of 
America, Vol. 67 (3), 823 - 826, 1980. 

9. F.J. Fahy, "Sound and Structural - Radiation Transmission 
and Response", Academic Press, 1985. 


18 




Material: Aluminium 
Density: 2710 Kg/m 3 
Elastic Modulus: 72 GN/M 
Thickness: 0.00635m 
Dimensions: a=1.0m, 

b=0.5m, 

Loss Factor: 0.01 


ure 1 . L— shaped plate structure showing plate 
characteristics. 
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Figire 3. Normalized power flow results for distributed loading. 
: Power input; — : power transfer; power ratio. 
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